Molecular Epidemiology and Diversity of SARS-CoV-2 in Ethiopia, 2020-2022 (preprint)

Ethiopia is the second most populous country in Africa and the sixth most affected by COVID-19 on the continent. Despite having experienced five infection waves, >499 000 cases, and ~7 500 COVID-19-related deaths as of January 2023, there is still no detailed genomic epidemiological report on the introduction and spread of SARS-CoV-2 in Ethiopia. In this study, we reconstructed and elucidated the COVID-19 epidemic dynamics. Specifically, we investigated the introduction, local transmission, ongoing evolution, and spread of SARS-CoV-2 during the first four infection waves using 353 high-quality near-whole genomes sampled in Ethiopia. Our results show that whereas viral introductions seeded the first wave, subsequent waves were seeded by local transmission. The B.1.480 lineage emerged in the first wave and notably remained in circulation even after the emergence of the Alpha variant. The B.1.480 was out-competed by the Delta variant. Notably, Ethiopia lack of local sequencing capacity was further limited by sporadic, uneven, and insufficient sampling that limited the incorporation of genomic epidemiology in the epidemic public health response in Ethiopia. These results highlight Ethiopia role in SARS-CoV-2 dissemination and the urgent need for balanced, near-real-time genomic sequencing.

Comparing SARS-CoV-2 sequencing methodologies during early phase detection of the Delta variant in South Africa (preprint)

Background: Genomic surveillance, with the aid of next-generation sequencing (NGS) technologies, revolutionized the SARS-CoV-2 pandemic. Coupled with high-performance analysis software, methodologies such as the Ion Torrent S5 and Illumina MiSeq dramatically improved the genomic surveillance capacity within South Africa during the height of the pandemic. Using de-identified remnant samples collected from Eastern Cape and analysis software, Genome Detective and NextClade, we compared the sequencing process, genomic coverage, quantification of mutations, and clade classification from sequence data generated by these two common “benchtop” NGS platforms. Results: Sequence data analysis revealed success rates of 175/183 (96%) and 172/183 (94%) on the Ion Torrent S5 and Illumina MiSeq, respectively. Internal quality metrics were assessed in terms of genomic coverage (>80%) and the number of mutations identified (<100). A greater number of higher-genomic coverage sequences were generated on the Ion Torrent S5 (99%) than on the Illumina MiSeq (80%) and <100 mutations was obtained by both platforms. Ion Torrent S5 generated high coverage sequences from samples having a broader range of viral loads (VL) compared to the Illumina MiSeq, which was less successful in sequencing samples with lower viral loads. Clade assignments were comparable across platforms which accurately differentiated between Beta (<45%) and Delta (≤30%) VOCs. A disparity in clade assignment was observed in <10% of sequences due to poor coverage obtained on the Illumina MiSeq, followed by a failure rate of ≤6% across the two platforms. Manual library preparation found both methods similar in terms of sample processing, handling of larger sample quantities, and clade assignment for SARS-CoV-2. Variability between the Ion Torrent S5 and Illumina MiSeq was observed in sequencing run duration (3,5 hrs vs 36 hrs), sequencing process (semi-automation vs manual), genomic coverage (99% vs 80%), and viral load requirements (broad range vs high VL). Conclusion: The Illumina MiSeq and Ion Torrent S5 are both reliable platforms capable of performing WGS with the use of amplicons and providing specific, accurate, and high throughput analysis of the SARS-CoV-2 whole viral genomes. Both sequencing platforms are feasible platforms for the genomic surveillance of SARS-CoV-2, each with its specific advantages and trade-offs.

Impact of intra-host immune adaptations on the evolution of SARS-CoV-2 S protein among individuals with SARS-CoV-2 infections in South Africa, 2020 to 2022 (preprint)

Background: Intra-host diversity studies are used to characterise mutational heterogeneity of SARS-CoV-2 infections to understand the impact of virus-host adaptations. This study investigated the frequency and diversity of the spike (S) protein mutations within SARS-CoV-2 infected South African individuals. Methods Single nucleotide polymorphism (SNP) assays and whole genome sequencing were performed on SARS-CoV-2 positive samples. Allele frequency (AF) was determined using TaqMan Genotyper software for SNP analysis and galaxy.eu for analysis of FASTQ reads. Results The SNP assays identified 5.3% (50/948) Delta cases with heterogeneity at delY144 (4%; 2/50), E484Q (6%; 3/50), N501Y (2%; 1/50) and P681H (88%; 44/50). Sequencing identified 9% (210/2381) cases with Beta, Delta, Omicron BA.1, BA.2.15, and BA.4 lineages with heterogeneity in the S protein. Heterogeneity was primarily identified at positions 19 (1.4%) with T19IR (AF 0.2-0.7), 371 (92.3%) with S371FP (AF 0.1-1.0), and 484 (1.9%) with E484AK (0.2-0.7), E484AQ (AF 0.4-0.5) and E484KQ (AF 0.1-0.4). Conclusion Mutations at heterozygous amino acid positions 19, 371 and 484 reduce recognition of neutralising antibodies, however the impact of the multiple substitutions at the same position is unknown. Therefore, we hypothesise that intra-host SARS-CoV-2 quasispecies with heterogeneity in the S protein facilitate competitive advantage of variants that can completely/partially evade host’s natural and vaccine-induced immune responses.

Efficacy of the AZD1222 (ChAdOx1 nCoV-19) COVID-19 Vaccine Against SARS-CoV-2 Variants of Concern (preprint)

In this South African phase 1/2b study, we demonstrated vaccine efficacy (VE) of two doses of AZD1222 for asymptomatic and symptomatic SARS-CoV-2 infection: 90.6% against wild-type and 77.1% against the Delta variant [≥]9 months after vaccination. VE against infection with the Beta variant, which preceded circulation of Delta, was 6.7%. Clinical trial identifierCT.gov NCT04444674

Continued Emergence and Evolution of Omicron in South Africa: New BA.4 and BA.5 lineages (preprint)

South Africa's fourth COVID-19 wave was driven predominantly by three lineages (BA.1, BA.2 and BA.3) of the SARS-CoV-2 Omicron variant of concern. We have now identified two new lineages, BA.4 and BA.5. The spike proteins of BA.4 and BA.5 are identical, and comparable to BA.2 except for the addition of 69-70del, L452R, F486V and the wild type amino acid at Q493. The 69-70 deletion in spike allows these lineages to be identified by the proxy marker of S-gene target failure with the TaqPath COVID-19 qPCR assay. BA.4 and BA.5 have rapidly replaced BA.2, reaching more than 50% of sequenced cases in South Africa from the first week of April 2022 onwards. Using a multinomial logistic regression model, we estimate growth advantages for BA.4 and BA.5 of 0.08 (95% CI: 0.07 - 0.09) and 0.12 (95% CI: 0.09 - 0.15) per day respectively over BA.2 in South Africa.

Omicron sub-lineages BA.4/BA.5 escape BA.1 infection elicited neutralizing immunity (preprint)

The SARS-CoV-2 Omicron (B.1.1.529) variant first emerged as the BA.1 sub-lineage, with extensive escape from neutralizing immunity elicited by previous infection with other variants, vaccines, or combinations of both. Two new sub-lineages, BA.4 and BA.5, are now emerging in South Africa with changes relative to BA.1, including L452R and F486V mutations in the spike receptor binding domain. We isolated live BA.4 and BA.5 viruses and tested them against neutralizing immunity elicited to BA.1 infection in participants who were Omicron/BA.1 infected but unvaccinated (n=24) and participants vaccinated with Pfizer BNT162b2 or Johnson and Johnson Ad26.CoV.2S with breakthrough Omicron/BA.1 infection (n=15). In unvaccinated individuals, FRNT50, the inverse of the dilution for 50% neutralization, declined from 275 for BA.1 to 36 for BA.4 and 37 for BA.5, a 7.6 and 7.5-fold drop, respectively. In vaccinated BA.1 breakthroughs, FRNT50 declined from 507 for BA.1 to 158 for BA.4 (3.2-fold) and 198 for BA.5 (2.6-fold). Absolute BA.4 and BA.5 neutralization levels were about 5-fold higher in this group versus unvaccinated BA.1 infected participants. The observed escape of BA.4 and BA.5 from BA.1 elicited immunity is more moderate than of BA.1 against previous immunity. However, the low absolute neutralization levels for BA.4 and BA.5, particularly in the unvaccinated group, are unlikely to protect well against symptomatic infection. This may indicate that, based on neutralization escape, BA.4 and BA.5 have potential to result in a new infection wave.

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance. (preprint)

Investment in Africa over the past year with regards to SARS-CoV-2 genotyping has led to a massive increase in the number of sequences, exceeding 100,000 genomes generated to track the pandemic on the continent. Our results show an increase in the number of African countries able to sequence within their own borders, coupled with a decrease in sequencing turnaround time. Findings from this genomic surveillance underscores the heterogeneous nature of the pandemic but we observe repeated dissemination of SARS-CoV-2 variants within the continent. Sustained investment for genomic surveillance in Africa is needed as the virus continues to evolve, particularly in the low vaccination landscape. These investments are very crucial for preparedness and response for future pathogen outbreaks.

Comparison of SARS-CoV-2 sequencing using the ONT GridION and the Illumina MiSeq (preprint)

Background: Over 4 million SARS-CoV-2 genomes have been sequenced globally in the past 2 years. This has been crucial in elucidating transmission chains within communities, the development of new diagnostic methods, vaccines, and antivirals. Although several sequencing technologies have been employed, Illumina and Oxford Nanopore remain the two most commonly used platforms. The sequence quality between these two platforms warrants a comparison of the genomes produced by the two technologies. Here, we compared the sequence quality produced by the Oxford Nanopore Technology GridION and the Illumina MiSeq for 28 sequencing runs. Results: Our results show that the MiSeq had a significantly higher number of sequences classified by Nextclade as good and mediocre compared to the GridION. The MiSeq also had a significantly higher sequence coverage and mutation counts than the GridION. Conclusion: Due to the low sequence coverage, high number of indels, and sensitivity to viral load noted with the GridION when compared to MiSeq, we can conclude that the MiSeq is more favourable for genomic surveillance, as successful genomic surveillance is dependent on high quality, near-whole genome sequences.

Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa (preprint)

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemic in southern Africa has been characterised by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, whilst the second and third waves were driven by the Beta and Delta variants respectively. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng Province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, predicted to influence antibody neutralization and spike function4. Here, we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

Management and containment of a SARS-CoV-2 Beta variant outbreak at the Malawi University of Science and Technology (preprint)

Outbreaks of COVID at university campuses can spread rapidly and threaten the broader community. We describe the management of an outbreak at a Malawian university in April 2021 during Malawi's second wave. Classes were suspended following detection of infections by routine testing and campus-wide PCR mass testing was conducted. Fifty seven cases were recorded, 55 among students, two among staff. Classes resumed 28 days after suspension following two weeks without cases. Just 6.3% of full-time staff and 87.4% of outsourced staff tested while 65% of students at the main campus and 74% at the extension campus were tested. Final year students had significantly higher positivity and lower testing coverage compared to freshmen. All viruses sequenced were beta variant and at least four separate virus introductions onto campus were observed. These findings are useful for development of campus outbreak responses and indicate the need to emphasize staff, males and senior students in testing.

Rapid replacement of the Beta variant by the Delta variant in South Africa (preprint)

The Beta variant of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in South Africa in late 2020 and rapidly became the dominant variant, causing over 95% of infections in the country during and after the second epidemic wave. Here we show rapid replacement of the Beta variant by the Delta variant, a highly transmissible variant of concern (VOC) that emerged in India and subsequently spread around the world. The Delta variant was imported to South Africa primarily from India, spread rapidly in large monophyletic clusters to all provinces, and became dominant within three months of introduction. This was associated with a resurgence in community transmission, leading to a third wave which was associated with a high number of deaths. We estimated a growth advantage for the Delta variant in South Africa of 0.089 (95% confidence interval [CI] 0.084-0.093) per day which corresponds to a transmission advantage of 46% (95% CI 44-48) compared to the Beta variant. These data provide additional support for the increased transmissibility of the Delta variant relative to other VOC and highlight how dynamic shifts in the distribution of variants contribute to the ongoing public health threat.

The continuous evolution of SARS-CoV-2 in South Africa: a new lineage with rapid accumulation of mutations of concern and global detection (preprint)

Global genomic surveillance of SARS-CoV-2 has identified variants associated with increased transmissibility, neutralization resistance and disease severity. Here we report the emergence of the PANGO lineage C.1.2, detected at low prevalence in South Africa and eleven other countries. The emergence of C.1.2, associated with a high substitution rate, includes changes within the spike protein that have been associated with increased transmissibility or reduced neutralization sensitivity in SARS-CoV-2 VOC/VOIs. Like Beta and Delta, C.1.2 shows significantly reduced neutralization sensitivity to plasma from vaccinees and individuals infected with the ancestral D614G virus. In contrast, convalescent donors infected with either Beta or Delta showed high plasma neutralization against C.1.2. These functional data suggest that vaccine efficacy against C.1.2 will be equivalent to Beta and Delta, and that prior infection with either Beta or Delta will likely offer protection against C.1.2.

Persistent SARS-CoV-2 infection and intra-host evolution in association with advanced HIV infection (preprint)

While most people effectively clear SARS-CoV-2, there are several reports of prolonged infection in immunosuppressed individuals. Here we present a case of prolonged infection of greater than 6 months with the shedding of high titter SARS-CoV-2 in an individual with advanced HIV and antiretroviral treatment failure. Through whole-genome sequencing at multiple time points, we demonstrate the early emergence of the E484K substitution associated with escape from neutralizing antibodies, followed by other escape mutations and the N501Y substitution found in most variants of concern. This provides support to the hypothesis of intra-host evolution as one mechanism for the emergence of SARS-CoV-2 variants with immune evasion properties.

A year of genomic surveillance reveals how the SARS-CoV-2 pandemic unfolded in Africa (preprint)

The progression of the SARS-CoV-2 pandemic in Africa has so far been heterogeneous and the full impact is not yet well understood. Here, we describe the genomic epidemiology using a dataset of 8746 genomes from 33 African countries and two overseas territories. We show that the epidemics in most countries were initiated by importations, predominantly from Europe, which diminished following the early introduction of international travel restrictions. As the pandemic progressed, ongoing transmission in many countries and increasing mobility led to the emergence and spread within the continent of many variants of concern and interest, such as B.1.351, B.1.525, A.23.1 and C.1.1. Although distorted by low sampling numbers and blind-spots, the findings highlight that Africa must not be left behind in the global pandemic response, otherwise it could become a breeding ground for new variants.

A novel variant of interest of SARS-CoV-2 with multiple spike mutations is identified from travel surveillance in Africa (preprint)

At the end of 2020, the Network for Genomic Surveillance in South Africa (NGS-SA) detected a SARS-CoV-2 variant of concern (VOC) in South Africa (501Y.V2 or PANGO lineage B.1.351)1. 501Y.V2 is associated with increased transmissibility and resistance to neutralizing antibodies elicited by natural infection and vaccination2,3. 501Y.V2 has since spread to over 50 countries around the world and has contributed to a significant resurgence of the epidemic in southern Africa. In order to rapidly characterize the spread of this and other emerging VOCs and variants of interest (VOIs), NGS-SA partnered with the Africa Centres for Disease Control and Prevention and the African Society of Laboratory Medicine through the Africa Pathogen Genomics Initiative to strengthen SARS-CoV-2 genomic surveillance across the region. Here, we report the first genomic surveillance results from Angola, which has had 21 500 reported cases and around 500 deaths from COVID-19 up to March 2021 (Supplemental Fig S1). On 15 January 2021, in response to the international spread of VOCs, the government instituted compulsory rapid antigen testing of all passengers arriving at the main international airport, in addition to the existing requirement to present a negative PCR test taken within 72 hours of travel. All individuals with a positive antigen test are isolated in a government facility for a minimum of 14 days and require two negative RT-PCR tests at least 48 hours apart for de-isolation, whilst all travelers with a negative test on arrival proceed to mandatory self-quarantine for 10 days followed by a repeat test. In March 2021, we received 118 nasopharyngeal swab samples collected between June 2020 and February 2021, a number of which were from incoming air travelers (Supplemental Fig S1). From these, we produced 73 high quality genomes (>80% coverage), 14 of which were known VOCs/VOIs (seven 501Y.V2/B.1.351, six B.1.1.7, one B.1.525), 44 of which were C.16 (a common lineage circulating in Portugal), and twelve of which were other lineages (Supplemental Fig S2). In addition, we detected a new VOI in three incoming travelers from Tanzania who were tested together at the airport in mid-February. The three genomes from these passengers were almost identical and presented highly divergent sequences within the A lineage (Figure 1A & 1B). The GISAID database contains nine other sequences reported to be sampled from cases involving travel from Tanzania, two of which are basal to the three sampled in Angola (Figure 1A, Supplemental Table S1). This new VOI, temporarily designated A.VOI.V2, has 31 amino acid substitutions (11 in spike) and three deletions (all in spike) (Figure 1C & 1D). The spike mutations include three substitutions in the receptor-binding domain (R346K, T478R and E484K); five substitutions and three deletions in the N-terminal domain, some of which are within the antigenic supersite (Y144{Delta}, R246M, SYL247-249{Delta} and W258L)4; and two substitutions adjacent to the S1/S2 cleavage site (H655Y and P681H). Several of these mutations are present in other VOCs/VOIs and are evolving under positive selection.

Escape of SARS-CoV-2 501Y.V2 variants from neutralization by convalescent plasma (preprint)

New SARS-CoV-2 variants with mutations in the spike glycoprotein have arisen independently at multiple locations and may have functional significance. The combination of mutations in the 501Y.V2 variant first detected in South Africa include the N501Y, K417N, and E484K mutations in the receptor binding domain (RBD) as well as mutations in the N-terminal domain (NTD). Here we address whether the 501Y.V2 variant could escape the neutralizing antibody response elicited by natural infection with earlier variants. We were the first to outgrow two variants of 501Y.V2 from South Africa, designated 501Y.V2.HV001 and 501Y.V2.HVdF002. We examined the neutralizing effect of convalescent plasma collected from six adults hospitalized with COVID-19 using a microneutralization assay with live (authentic) virus. Whole genome sequencing of the infecting virus of the plasma donors confirmed the absence of the spike mutations which characterize 501Y.V2. We infected with 501Y.V2.HV001 and 501Y.V2.HVdF002 and compared plasma neutralization to first wave virus which contained the D614G mutation but no RBD or NTD mutations. We observed that neutralization of the 501Y.V2 variants was strongly attenuated, with IC50 6 to 200-fold higher relative to first wave virus. The degree of attenuation varied between participants and included a knockout of neutralization activity. This observation indicates that 501Y.V2 may escape the neutralizing antibody response elicited by prior natural infection. It raises a concern of potential reduced protection against re-infection and by vaccines designed to target the spike protein of earlier SARS-CoV-2 variants.

Emergence and rapid spread of a new severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) lineage with multiple spike mutations in South Africa (preprint)

Continued uncontrolled transmission of the severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) in many parts of the world is creating the conditions for significant virus evolution. Here, we describe a new SARS-CoV-2 lineage (501Y.V2) characterised by eight lineage-defining mutations in the spike protein, including three at important residues in the receptor-binding domain (K417N, E484K and N501Y) that may have functional significance. This lineage emerged in South Africa after the first epidemic wave in a severely affected metropolitan area, Nelson Mandela Bay, located on the coast of the Eastern Cape Province. This lineage spread rapidly, becoming within weeks the dominant lineage in the Eastern Cape and Western Cape Provinces. Whilst the full significance of the mutations is yet to be determined, the genomic data, showing the rapid displacement of other lineages, suggest that this lineage may be associated with increased transmissibility.

High Resolution analysis of Transmission Dynamics of Sars-Cov-2 in Two Major Hospital Outbreaks in South Africa Leveraging Intrahost Diversity (preprint)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes acute, highly transmissible respiratory infection in both humans and wide range of animal species. Its rapid spread globally and devasting effects have resulted into a major public health emergency prompting the need for methodological interventions to understand and control its spread. In particular, The ability to effectively retrace its transmission pathways in outbreaks remains a major challenge. This is further exacerbated by our limited understanding of its underlying evolutionary mechanism. Using NGS whole-genome data, we determined whether inter- and intra-host diversity coupled with bottleneck analysis can retrace the pathway of viral transmission in two epidemiologically well characterised nosocomial outbreaks in healthcare settings supported by phylogenetic analysis. Additionally, we assessed the mutational landscape, selection pressure and diversity of the identified variants. Our findings showed evidence of intrahost variant transmission and evolution of SARS-CoV-2 after infection These observations were consistent with the results from the bottleneck analysis suggesting that certain intrahost variants in this study could have been transmitted to recipients. In both outbreaks, we observed iSNVs and SNVs shared by putative source-recipients pairs. Majority of the observed iSNVs were positioned in the S and ORF1ab region. AG, CT and TC nucleotide changes were enriched across SARS-COV-2 genome. Moreover, SARS-COV-2 genome had limited diversity in some loci while being highly conserved in others. Overall, Our findings show that the synergistic effect of combining withinhost diversity and bottleneck estimations greatly enhances resolution of transmission events in Sars-Cov-2 outbreaks. They also provide insight into the genome diversity suggesting purifying selection may be involved in the transmission. Together these results will help in developing strategies to elucidate transmission events and curtail the spread of Sars-Cov-2

Major new lineages of SARS-CoV-2 emerge and spread in South Africa during lockdown. (preprint)

In March 2020, the first cases of COVID-19 were reported in South Africa. The epidemic spread very fast despite an early and extreme lockdown and infected over 600,000 people, by far the highest number of infections in an African country. To rapidly understand the spread of SARS-CoV-2 in South Africa, we formed the Network for Genomics Surveillance in South Africa (NGS-SA). Here, we analyze 1,365 high quality whole genomes and identify 16 new lineages of SARS-CoV-2. Most of these unique lineages have mutations that are found hardly anywhere else in the world. We also show that three lineages spread widely in South Africa and contributed to ~42% of all of the infections in the country. This included the first identified C lineage of SARS-CoV-2, C.1, which has 16 mutations as compared with the original Wuhan sequence. C.1 was the most geographically widespread lineage in South Africa, causing infections in multiple provinces and in all of the eleven districts in KwaZulu-Natal (KZN), the most sampled province. Interestingly, the first South-African specific lineage, B.1.106, which was identified in April 2020, became extinct after nosocomial outbreaks were controlled. Our findings show that genomic surveillance can be implemented on a large scale in Africa to identify and control the spread of SARS-CoV-2.